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We investigate the fermionic SU(A^) Hubbard model on the two-dimensional square lattice for weak 
to moderate interaction strengths using one-loop renormalization group and mean-field methods. 
For the repulsive case (7 > at half filling and small A'^ the dominant tendency is towards breaking 
of the SU(A'^) symmetry. For > 6 staggered flux order takes over as the dominant instability, 
in agreement with the large-A'' limit. Away from half filling for A' = 3 the system rearranges the 
particle densities such that two flavors remain half filled by cannibalizing the third fiavor. In the 
attractive case and odd A' a full Fermi surface coexists with a superconductor in the ground state. 
These results may be relevant to future experiments with cold fermionic atoms in optical lattices. 



Introduction: After the celebrated observation of 
Bose-Einstein condensation ultracold atom systems 
receive growing attention in the field of condensed mat- 
ter physics. Recently, also quantum degenerate Fermi 
gases have been realized 0i S S 13 j opening up the pos- 
sibility to study phenomena such as BCS superfluidity 
in a new context. As a further important advance, op- 
tical lattices have been used to realize the transition be- 
tween a bosonic superfluid and a Mott insulator P|. It 
has thus been demonstrated that cold atoms systems can 
become a very flexible and clean laboratory for many ex- 
citing phenomena from the purview of condensed mat- 
ter or interacting many particle systems. In particular, 
it has been suggested Q that cold fermions in optical 
lattices may help to understand the notorious complexi- 
ties of strongly correlated solid state systems such as the 
cuprate high-temperature superconductors. 

Besides the realization of phenomena that are known 
to exist in some form in condensed matter systems, it 
is also interesting to ask whether the degrees of freedom 
offered by cold atoms could give rise to states of mat- 
ter that do not have obvious counterparts in the physics 
of interacting electrons. Typical electron systems, at 
least in the flrst approximation, possess SU(2) spin ro- 
tational symmetry which can be broken spontaneously, 
leading to magnetic phenomena such as ferro- and anti- 
ferromagnetism. For alkali atoms, the nuclear spin / and 
electron spin S are combined in a hyperfine state. Its 
total angular momentum F can be different from 1/2, 
and for each F there are 2F + 1 hyperfine states differ- 
ing by their azimuthal quantum number mp- E.g. for 
the fermionic '^''K, the nuclear spin is / = 4 and the low- 
est hyperfine multiplet (at weak fields) has F = 9/2. In 
magnetic traps only a subset of these 2F + 1 states (the 
low-field-seekers) can be trapped j^, but this constraint 
can be avoided by using all-optical traps . 

In fact, coexistence of the three hyperfine states \F = 
9/2, mp = -5/2,-7/2,-9/2) of in an optical trap 
has already been realized, with tunable interactions due 
to Feshbach resonances between itif = —5/2/ — 9/2 and 
mp = —7/2/ — 9/2, respectively A situation with 
strong attractive interaction between all three compo- 
nents can be realized e.g. for the spin polarized states 



with rris = 1/2 in ^Li where the tri ple t scattering length 
a = — 2160ao is anomalously large [lOl |. 

Optical lattices are created by a standing light wave 
leading to a periodic potential for the atomic motion of 
the form V{x) = Vq J^i cos^ {kxi) where k is the wavevec- 
tor of the laser, i labels the spatial coordinates and the 
lattice depth Vq is usually measured in units of the atomic 
recoil energy En = fP'k^ /2m. In the following we will 
consider the 2D case where i — \,2. It has been shown 
11] that the Hubbard model with a local density-density 
interaction provides an excellent description of the low- 
energy physics. Here we are interested in a situation 
where fermionic atoms with N different spin states ( "fla- 
vors") m are loaded into the optical lattice. We thus 
consider a Hubbard Hamiltonian 
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Here = J^i,™ is the total number density of atoms 
on site i which can be written in terms of creation and 
annihilation operators according to nj.„j — (^imCi^m- The 
interaction (second term in Eq. ^ is invariant under lo- 
cal U(A^) rotations of the N flavors with different m. The 
hopping term of the atoms between nearest neighbors {ij) 
reduces the invariance of the complete Hamiltonian to a 
global U(Af) symmetry. Stripping off the overall U(l) 
phase factor, we arrive at the SU(A^) Hubbard model. 
Note that in the optical lattice the effective Hubbard pa- 
rameters are given hy t — En (2/ ^Jii) ^ exp (—2^^) and 

\J = Enask^J^pK^ where ^ = {Vq/Ee)^^^ and is the 
s-wave atomic scattering length. 

The fermionic SU(iV) Hubbard model for U > on 
the two-dimensional (2D) square lattice was studied in 
the large-A^ limit tl.2] in the early days of high-T^ su- 
perconductivity, mainly as a controllable limit connected 
to the then physically relevant case N = 2. A gener- 
alized SU(Af) model could describe orbitally degenerate 
electronic states in crystals, but it is likely that in these 
systems different overlaps between the orbitals pointing 
in distinct lattice directions break the SU(A^) invariance. 

In the following we focus on the density region near 
half band flUing with an average A^/2 fermions per site. 



In the conventional N = 2 Hubbard model at half filling 
the ground state exhibits spin-density wave (SDW) order, 
and when the filling is changed d-wave superconductiv- 
ity is very likely '13'|. The SDW state breaks the trans- 
lational invariance and spin-up and spin-down electrons 
(for staggered moment along the z-dircction) occupy the 
two sublattices differently (see Fig. For large N and 
small exchange interactions J, staggered flux order is ex- 
pected to dominate over the SU(7V)-breaking states [T^. 

One-loop renormalization group for the half-filled band 
and general N and U > 0: First let us analyze the one- 
loop renormalization group (RG) flow for the half filled 
band. We apply the perturbative temperature-fiow RG 
method of Ref. 14] that has proved to give good results 
for N — 2. As initial condition one fixes the interaction 
at a high temperature of the order of the bandwidth. 
Then the RG fiow describes the change of the interactions 
as the temperature is lowered and perturbative correc- 
tions due to one-loop particle-hole and particle-particle 
processes are taken into account. The interaction is de- 
scribed by a coupling function V^(fci,fc2,^3) |l5j| . where 
the fiavor indices mi and belonging to the first in- 
coming particle with wavevector ki and the first outgo- 
ing particle with are the same. Similarly m2 — m^. 
The second outgoing wavevector k^ is fixed by momen- 
tum conservation on the lattice. As in the N ^ 2 case 
the RG flow goes to strong coupling. This means, as 
we start the flow at high temperatures with a purely lo- 
cal interaction V{ki,k2,k3) = [/, some coupling func- 
tions start to grow when the temperature is reduced and 
flnally leave the perturbative range. At this tempera- 
ture scale we stop the RG flow and analyze which class 
of coupling constants grows most strongly towards low 
T. In close analogy with the spin- 1/2 case we consider 
couplings in the charge channel Vc{k, k', q) = NV{k + 
q,k',k) — V{k\k -\- q,k) and in the SU(iV) symmetry 
breaking channel, Vs{k, k', q) = —V{k', k-\-q, k), which, if 
divergent, signals a singular response for a small external 
fleld coupling to one of the iV^ — 1 generators of SU(A^). 
We deflne averages over the Fermi surface, V^^^{Q) — 

§pgd(t)k §ps, d(pk' 9t{k)9e{k')Vc/s{k,k',q)with gi{k) = 1 
in the s-wave channel or .92 (fc) = (cosfc^, — cosfc,y)/\/2 in 
the d-wave channel. For N ^ 2 the couplings Vs{Q) in the 
s-wave channel with Q = (vr, tt) diverge most strongly, 
indicating the formation of an antiferromagnetic (AF) 
SDW state. Note that for larger N the charge couplings 
can become quite large as they scale with N. 

Analyzing the RG flows to strong coupling, the picture 
is as follows: for N < 6 the s-wave coupling function 
V^{Q) is the strongest divergent family of coupling con- 
stants, signalling a dominant tendency towards breaking 
of the SU(7V) symmetry with staggered two-sublattice 
real space dependence. For N > 2 (especially for odd 
N) this leads to the interesting question how the N/2 
particles per site will arrange themselves on the bipar- 
tite square lattice (see Fig. Below we describe what 
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FIG. 1: a); AF spin-density wave state for TV = 2. Spin- 
up and spin-down particles occupy the two sublattices with 
different probabilities (here idealized to and 1). b) SU(3)- 
breaking flavor-density wave state for N — 3. Flavors 1 and 
2 prefer one sublattice, flavor 3 the other, c) Staggered flux 
state for A*' > 6: the arrows indicate the particle currents. 
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FIG. 2: Left plots; RG results for SU(3) at half filling and 
U = 4t. Upper left; Flow of the coupling constants Vg" in the 
SU(3)-breaking channel (dashed line) and Vc in the staggered 
flux (SF) channel, averaged around the Fermi surface (FS). 
Lower left; V3{k,k' ,Q) (dashed line) and Vc{k,k' ,Q) with k 
fixed on the FS at (tt, 0) and k' moving around the FS after 
termination of the flow. Right plots; The same for SU(8) at 
half filling. 



happens in a mean-field analysis. 

For TV > 6 the flow to strong coupling changes qual- 
itatively. Now the leading divergence is in the charge 
couplings Vc(fc, k', Q) with a d^-^-y-^-wave dependence on 

k and k' . V/J diverges more strongly than (see Fig. |21l, 
albeit at lower temperature T « 0.014i for U = At and 
N = 7. This signals a tendency towards staggered flux 
(SF) order with long range ordering of the expectation 
value $SF = Efc,™(cosA::, - cos fcj,)(ct ^c^^^ „J. This 
result agrees with the large-TV limit for small exchange 
interactions J [l3|- The SF state has surfaced several 
times for the SU(2) case in connection with the high- Tc 
cuprates and related models (1^ il6||J, also as d-density 
wave state (although the particle density is not modu- 
lated). Its quasiparticles have a wavevector-dependent 
energy gap that vanishes at A: = (±7r/2, ±7r/2). Nonzero 
^SF breaks translational and time-reversal symmetry 
with alternating particle currents around the plaqucttes 
(see Fig. If the particles were charged, their motion 
would give rise to alternating magnetic moments point- 
ing out of the plane, hence the name staggered flux state. 
Note that ^sf is SU(TV)-invariant and no continuous 
symmetries are broken. Correspondingly the SF state 
can order at finite temperatures in 2D. For the same rea- 
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son it may be possible that the staggered flux state sets 
in for somewhat lower N than the critical A'^ = 6 in our 
one-loop RG study that neglects collective fluctuations. 

Away from half filling the flow of the dominant (tt, tt)- 
instability gets cut off at some low energy scale that in- 
creases with the distance to half filling. Below that scale 
there is a tendency towards da;2_j,2-wave Cooper pair- 
ing. However the energy scales for pairing instabilities 
become very small with increasing N. Below we discuss 
superfluid pairing for iV > 2 in the attractive case U < 0. 

Ground state near half filling for N = 3: Having es- 
tablished that SU(A'^) symmetry breaking at wave vector 
(tt, tt) is the dominant instability of the repulsive model 
near half filling with iV < 6, we now turn to a mean-field 
description of the ground state for = 3. We decouple 
the interaction terms in the particle-hole channel with 
local mean-fields (cj^, ^cp^i) = M^fj^i. The hermitian local 
mean-field matrix Map can be decomposed into a trace- 
ful part Mq proportional to the identity matrix to and a 
traceless part X]a=i 8 ^'^""^a with the 8 generators ta of 
the fundamental representation of SU(3). A finite value 
of one of the traceless components breaks the SU(3) in- 
variance. We will now restrict the analysis to commen- 
surate order, where only uniform and staggered compo- 
nents of a commuting subset of the 9 M'^ acquire nonzero 
expectation values. SU(3) has rank 2 and the two diago- 
nal generators commute mutually and with the identity 
matrix. We can choose these three degrees of freedom to 
be contained in the three flavor density mean-fields {ua). 

The results of T = mean-field solutions are shown in 
Fig. 01 At half filhng, n = 1.5/site, the SU(3)-breaking 
creates a flavor-density wave: two flavors prefer one sub- 
lattice with equal density, while the third flavor goes 
predominantly on the other sublattice with a somewhat 
larger density modulation. The staggered components 
do not add up to zero. Thus there is a charge density 
wave accompanying the SU(3) symmetry breaking. For 
[/ = 3i the mean-field Tc for this state is ~ 0.45t, but in 
the one-loop RG is it is reduced down to ~ 0.12i. 

Fig. 121 also describes the results away from half filling. 
For example at J7 = l.&t and 1.42 < n < 1.48/site, two 
flavors order with opposite staggered densities on the two 
sublattices, keeping their individual average density at 
half filling. Since the total density is less than half filling, 
the third fiavor gets decimated with uniform density of 
(n — 1). As can be seen from the right plot in Fig. (21 this 
state only occurs above a critical interaction strength Uc 
that increases from zero with increasing distance to rt = 
1.5/site. Note that the depletion of one flavor allows the 
system to preserve the commensurate order away from 
commensurate band filling. A similar pinning of a part 
of the system to half filling is found in ladder systems 
01 . We add that for larger U and close to half filling 
(1.46 < n < 1.5/site for U = At) we find a another 
regime where the mean field equations converge slowly 
and microscopic phase separation might occur. 

Attractive SU(N) Hubbard model: Next let us con- 
sider the attractive interactions U < 0. In the SU(2) 
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FIG. 3: Left: Uniform densities of the = 3 flavors (solid 
line: flavor 1 and 2, crosses: flavor 3) and staggered densities 
(squares) vs. total density for U = 1.6t. Right: Difference 
in uniform density between the two majority flavors and the 
minority flavor versus interaction U and total density (n) per 
site. The scalebar indicates the density difference per site. 



case in 2D, there is a power-law s-wave singlet super- 
conductor/superfluid (SSC) below a Kosterlitz-Thouless 
transition away from half filling At half filling 

the Kosterlitz-Thouless Tc is zero, and SSC and charge- 
density wave (CDW) mean field states are degenerate 
and in the true ground state, both types of order coex- 
ist. One-loop RG finds in this case that CDW and SSC 
susceptibilities are perfectly degenerate and diverge to- 
gether at low T. This symmetry is destroyed for larger 
N > 2 and the CDW susceptibility grows much faster 
than the one for SSC. Correspondingly we expect the 
ground state to have CDW long range order only. This 
is corroborated by a mean-field calculation for the SU(3) 
case that shows that the CDW order suppresses any kind 
of SSC admixture, and the CDW ground state energy is 
lower than that of the SSC state. 

We now consider the generic case sufficiently far away 
from half filling. Then the dominant instability is on- 
site pairing. We decouple the interaction as Hu^^f = 
— h fl ct - A«ct + h.c. with the local mean-fields 
Aa/3 = -U I]fc(CfeQC_j^) = -A/3a. For TV > 2 these even 
parity gap functions Aq,/? transform non-trivially under 
SU(A^). Depending on the global gauge, Aq,^ takes dif- 
ferent values. This is unlike the SU(2) case where even 
parity gap functions are singlets and invariant under spin 
rotations [l^- For SU(2) the ground state is degenerate 
with respect to the global phase of the gap function, and 
long- wavelength variations of the latter are gapless in ab- 
sence of long-range forces. In the SU(Af) case we find 
a higher degeneracy and more gapless modes. It turns 
out that for SU(3) all gap functions with with the same 
Aq = X]a/3 ^'^^ degenerate and have the same to- 

tal density of states. Apart from the global phase there 
are four additional gapless modes, two associated with 
the internal phases between A12, A13 and A23, and two 
modes modulating IA12I, IA13I and IA23I with fixed Aq. 

A particularly simple choice in the degenerate man- 
ifold is A12 = Aq and A13 — A23 = 0. Then flavor 
3 remains completely unpaired and metallic. Since we 
can always rotate into this gauge, all SU(3) s-wave su- 
perconducting mean-field states are one-third (neutral) 
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metals and two-thirds superfluids. The gauge with only 
Ai2 ^ makes the symmetry breaking pattern obvious. 
The original symmetry group of the problem SU(3) (E) 
U(l) with nine generators gets broken down to an SU(2) 
in flavor 1 and 2, leaving A12 invariant, and an addi- 
tional U(l) that acts on the phase of the unpaired flavor 
3. This leaves 5 generators broken, yielding the collective 
modes described above. For 3/8 band filling and U = 4t, 
the mean-field Tc is ~ O.lTt. The coexistence of a full 
Fermi surface with a superconductor should have inter- 
esting consequences. For example the collective modes 
may be subject to damping below twice the gap frequency 
and could hence render the ungapped fermionic spectrum 
observable. Experimentally, these Goldstone modes can 
be detected in the spectrum of elementary excitations 
measured via Bragg scattering off two non-coUinear laser 
beams with frequency and momentum difference w,g |2C| . 
By monitoring the number of scattered atoms as a func- 
tion of uj, this technique yields the dynamical structure 
factor S{q,uj). The collective modes will then lead to 
peaks in the scattering cross section. 

Theoretically, an additional weak p-wave attraction 
could trigger a superfluid transition of the unpaired fla- 
vor at much lower temperatures, leading to a coexistence 
of even- and odd-parity superfluidity. 

The SU(4) case is more complicated. There the degen- 
eracy of the ground state is subject to more constraints 
than just constant Aq. The mean- field solutions have 
IA12I = IA34I, IA13I = IA24I and |Ai4| = IA23I. The 
single particle spectrum is fully gapped. 

Conclusions: The fermionic SU(A^) Hubbard model 
on the 2D square lattice can possibly be realized with 
ultracold atoms in an optical lattice. Its ground states 
may exhibit phenomena that do not occur right away in 
traditional solid state systems. For U > we find a stag- 
gered flux state for large > 6 at half band fllling where 
the particles run around the plaquettes of the lattice in 



an alternating way. This state has a partially gapped 
excitation spectrum with nodes along the Brillouin zone 
diagonals, which may be detectable via the momentum 
distribution function. Near half fllling for iV = 3 we find 
a redistribution of the particle densities where two of the 
three flavors remain half filled and occupy different sub- 
lattices while the third flavor becomes depleted. Finally, 
in the attractive case U < we point out that the s- 
wave paired superfluid states may exhibit new collective 
modes. For A'^ = 3 a third of the particles remains un- 
gapped, leading to a full Fermi surface coexisting with 
the superfluid. We expect this to be a general feature for 
odd N, also in three dimensions or in absence of a lattice 
potential. 

Finally, a comment on the temperature scales which 
we have given in terms of the hopping parameter t. It 
has been shown Q that if the optical lattice is switched 
on slowly after termination of evaporative cooling, an 
additional adiabatic cooling process takes place. The fi- 
nal temperature is given by the identity Tinitiai/TV^free ~ 
Tfinai/JFjattice where Tp,free(iattice) denote the Fermi tem- 
perature of the free atomic cloud and in the presence of 
the lattice, respectively. In particular, in 2D at half filling 
one has Tpiattice = 4t. As a result, the critical atomic 
temperatures which have to be reached before the lat- 
tice is switched on can be obtained from our results via 
the substitution t Tp^hcc/'^- For the s-wave superfluid 
phase [U < 0) and the flavor-density wave states {U > 0) 
we therefore find transition temperatures of order Q.OBTp 
which are within experimental reach. 
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